Higher order antibunching is not a rare phenomenon 
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Abstract 

Since the introduction of higher order nonclassical effects, higher order squeezing has been reported in a 
number of different physical systems but higher order antibunching is predicted only in three particular cases. In 
the present work, we have shown that the higher order antibunching is not a rare phenomenon rather it can be 
seen in many simple optical processes. To establish our claim, we have shown it in six wave mixing process, four 
' wave mixing process and in second harmonic generation process. 
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1 Introduction 

^ Antibunching and squeezing do not have any classical analogue Higher order extensions of these nonclassical states 

have been introduced in recent past Among these higher order nonclassical effects, higher order squeezing has already 

been studied in detail [??,|S1IH1E| but the higher order antibunching (HOA) is not yet studied rigourously. The idea of 
HOA was introduced by Lee in a pioneering paper [O] in 1990, since then it is predicted in two photon coherent state [S], trio 
y—i coherent state JOI and in the interaction of intense laser beam with an inversion symmetric third order nonlinear medium 
J>" Jl]. From the fact that in last 15 years HOA is reported only in three particular cases, HOA appears to be a very rare 
phenomenon. The present study aims to establish that this apparent rarity is not due to any physical reason. To establish 
that, we have shown the existence of HOA in six wave mixing process, four wave mixing process and in second harmonic 
QQ generation process. 

Using the negativity of P function Lee introduced the criterion for HOA as 

in ■ 

^ where N is the usual number operator, (TV^*^) = {N{N — 1)...{N — i + 1)) is the ith factorial moment of number operator, 
{) denotes the quantum average, I and m are integers satisfying the conditions I < m < 1 and the subscript x denotes a 
^^articular mode. Ba An QH] choose m — 1 and reduced the criterion of Zth order antibunching to 

X 



A.,i = ' - K (2) 



Nt''^)<(N^'A{N^). (3) 



Physically, a state which is antibunched in Ith order has to be antibunched in (/ — l)th order. Therefore, we can further 
simplify (pj as 

(7Vi'+i)) < (7V«) (N.) < (iVi'-i)) {N^f < {Nf') {N^f < ... < {N^r' (4) 
and obtain the condition for I — th order antibunching as 

d(0 = (iVi'+^))-(iV.)'+'<0. (5) 

This simplified criterion ^ coincides exactly with the physical criterion of HOA introduced by Pathak and Garica [llj . 
Here we can note that d(l) = and d(l) > corresponds to higher order coherence and higher order bunching (many 
photon bunching) respectively. As we have already mentioned, higher order antibunching is not yet studied rigourously and 
apparently, higher order antibunching is very rare since it is reported only in three particular cases in last 15 years. The 
present work aims to show that its not really a rare phenomenon rather it can be seen in many simple optical processes. To 
establish that we have used the criterion Q and short time approximated solutions of equation of motions corresponding 
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to various Hamiltonians (such as six wave mixing process, four wave mixing process and second harmonic generation) and 
have shown that HOA can be seen in all the physical systems selected for the present study. In the next section we present 
a second order operator solution of the equation of motion of six wave mixing process in detail and use that to show the 
existence of HOA in six wave mixing process. In section 3 and 4 we have studied the possibilities of observing higher 
order antibunching in four wave mixing process and in second harmonic generation process respectively. Finally section 5 is 
dedicated to conclusions. 



2 Six wave mixing process 

Six wave mixing may happen in different ways. One way is that two photon of frequency tui are absorbed (as pump photon) 
and three photon of frequency UJ2 and another of frequency ui^ are emitted. The Hamiltonian representing this particular six 
wave mixing process is 

H = a^auji + b'^buj2 + c^ccja + g{a)'^b^c + a^h'^'^c^) (6) 

where a and a^are creation and annihilation operators in pump mode which satisfy [a, a^^]=l, similarly 6, 6^ and c, are 
creation and annihilation operators in stokes mode and signal mode respectively and g is the coupling constant. Substituting 
A = ae*"i*, B = 6e*'^^* and C = ce*"^* ^e can write the Hamiltonian ^ as 

H = A^Auji + B^Buo2 + C^Cuo3 + g{A^^B^C + A^B^^C^) (7) 

Since the Hamiltonian is known, we can use Heisenberg's equation of motion (with ft = 1): 

^ = ^ + *[^,^] (8) 
and short time approximation to find out the time evolution of the essential operators. From equation Q we have 

[H,A] = -Auoi-2gA'^B^C. (9) 



From ^ and ^ we have 
Similarly 
and 

C = -igA^B^"^. (12) 
We can find the second order differential of A using (jHJ and Hlflll2|l as 



A = iAuji - iAui - i2gA^B^C = -2igA^B^C. (10) 
13 = -iigA^B'^^C'^ (11) 



A = ^ + i[H,A] 

= Ag'^AB^'^B^C^'C - \%g^A'^A^B^'^B^0C - i&g'^A'^A^B^BC^C 

- 2g^A^A^B^^B^ - ISg^A^ A^B^^B^ ~ 36g^A^A^B^B 

- Ug^A'iA^C^C - 12g2AtA2. 



(13) 



Now by using the Taylor's series expansion 
and substituting lfT?1|l and (tT3|l in (tlHl we get 



(14) 



A{t) = A~ 2igtA'^B^C 

+ gH^ [2AB^'^B^C^C - gA^i A^B^^B'^C^'C - ISAtA^fit^c'tc (15) 
- A^A^B'i^B^ - 9A^A^B^^B^ - ISA'iA^B^B - GA^A^C^C ~ 6A^A^] 

or, 

A{t) = A~ 2igtA^B^C 

+ gH^ [ZAB'i^B^Nc - 9NaAB^^B^Nc - ISNaANbNc (16) 
- NaAB^^B^ - QNaAB^^B^ - ISNaANb - GNaANc - 6NaA] 

where Na = A'' A, Nb = B^B, Nc = C^C^. The Taylor series is valid when t is small, so this solution is valid for a short 
time and that is why it is called short time approximation. The above calculation is shown in detail as an example. Following 

■'A short time approximated expression of time evolution of annihilation operator in pump mode of six wave mixing process described by J6j is 
also derived in [2] but unfortunately their solution contain some mistakes. 
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the same prescription, we can find out the time evolution of B and C or any other creation and annihilation operator that 
appears in the Hamiltonian of matter field interaction. This is a very strong technique since this straight forward prescription 
is vaHd for any optical process where interaction time is short. After obtaining the analytic expression for time evolution of 
annihilation operator, now we can use it to check whether it satisfies condition or not. 

Let us start with the possibility of observing first order antibunching. From equation (tT5)l , we can derive expressions for 
N{t) and N^^^t) as 

N{t) = A^A-2tgt{A^^B^C - A^Bl^C'<) 

+ gH^ [8A'<AB'<^B^C'<C - 18A''^A'^B'<^B^C'<C - 36A'<^A^B^BC'<C - 2A^^A^B^^B^ (17) 
- 18At2A2Bt2s2 _ 36At2A2BtS - UA^^A^C^C + AB^^B^C^C - UA^^A^] 

and 

iV(2)(i) = At2(t)A2(t) = ^t2^2 _ 2igt (2At3AS3C - A^^B^'C - 2AtA3st3Ct - A^B^^C^) 
+ gH"^ [2AA''^A^B'i^B^C'fC + 32A^ AB^f^B^C^C + AB^^B^C^C 

- SGA^^A^B^^B^C^fC - 72A^^A^B'<BC'^C - ISA^^A^B^^B^C^C 

- 36At2A2^t5C'tC-4.4t3^35t353_3g^t3^35t2B2 (18) 

- 72At3^3stB - 2At2A2Bt3B3 - ISA^^A^B^^B^ 

- 4At4B*^C2 - AA*B^^C^^ - 36At2^2^ts 

- 24A^^A^C^C - UA^^A^C'fC - 24A^^A^ - 12At2A2] . 

In the present study, we have taken all the expectations with respect to |q; > |0 > |0 > for simplification. This assumption 
physically means that initially a coherent state (say, a laser) is used as pump and before the interaction of the pump with 
atom, there was no photon in signal mode (6) or stokes mode (c). Thus the pump interacts with the atom and causes 
excitation followed by emissions. Now from (I17|l and Ijl8|l . we have 

{N{t)f = \a\''~24gV\a\^. (19) 



{N\t)) = - .g2t2 (24|a|6 + 12\a\^) (20) 

where A\a >~ a\a >. Now by using Ijl9|l and Ij20|l we can show that the six wave mixing process satisfies the criterion of 
antibunching © since: 

d(l) = {N(^Ht))-{Nit)f 

= [|a|4 + gH^{-24\a\^ - 12\a\*)] - [\a\^ - 2452t2|„|6] (21) 

= -1252i2|^|4, 

From the last equation it is clear that d(l) is always negative, i.e. it always shows usual antibunching. Essentially, this is 
a nonclassical state but mere satisfaction of nonclassicality or antibunching is not enough because we are looking for HOA. 
Let us see what happens in the next higher order. 

For the study of second order of antibunching, A'^{t) can be obtained by using and operator ordering techniques: 

^3 (t) = A3 - 2igt {SA^A^B'^C + 3AB^C) 

+ gH^ [GA^B^^B^C^C - 27A^A^B^^B^C^C - 5AA'<A'^B^BC^C 

- 54A3stBCtC - 27A3Bt252(;;tc - 3AtA4Bt3B3 _ 27 A-^A^B^^B^ - 5AA'<A'^B'<B (22) 

- 3A3Bt353 _ 27yl3Bt2B2 - ISvlU^CtC - ISA^C^C 

- 54A^B'<B - 12At2AB6c2 - UA^B^C^ - ISA^A^ - 18A^] 

Then A^^{t) can simply be written as, 

At3(t) = At3 + 2^5^ (3At2ASt3ct + SAtSts^t) 

+ ^2^2 [6At3st3s3CtC - 27v4t4^_Bt252(^tc - bAA'f^ AB^ BC^ C 

- 54y4t3st5ctC - 27At35t252(^tc - 3^+4^5+3^3 _ 27^+4^5+252 _ 54^14^515 (23) 

- 3^+35+353 - 27.4+35+252 - 18y4t'*AC+C - 18A+3C+C 

- 54.4+35+5 - 12y4+y425+6C+2 - 12.45+6(7+2 - ISA+^A - 18A+3] 

Last two equations can be used to calculate the third factorial moment (N^'^^t)) of the number operator N as 
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7V(3) (t) = At3^3 - 2igt {SA^^A^B^C + SA^^AB^C - SA^^A^B^'^C^ - SA^A^B^'^C^) 
+ gH^ [48A^^A^B^^B^C^C - MA^^A^B^^B'^C^C - bAA^^A^B^^B'^C^C 

- 108At3A3BtBCtC+108At2A2Bt3B3CtC + 36yltABt3s3c'tc 

- 6At4A4Bt3B3 _ 54At3A3Bt2s2 - 108At3A3stS - 12At5AB6c2 - 12AtA5Bt6ct2 

- 12At4B6c2 - 12A4Bt6ct2 - 36At3A3CtC - SGAt^A^CtC - SGAt^A^ - 36At3A3] . 

Taking the expectation value with respect to the intial state we obtain 

(Ar(3)(t)) ^ |a|6 -g2^2 ^36|a|8 + 36|a|6) . (25) 

On the other hand, 

{N{t)f = \a\'^ -mgh''\a\^ (26) 

By using the last two equations we can see that the pump mode of six wave mixing process satisfy the criteria of antibunching 
of second order @. Since, 

d(2) = (iV(3)(t))-(7V(t))3 

= - g^e{m\a\^ + 36|a|6)] - - mgH''\a\^] (27) 

= -3652i2|a|6 

is always negative. 



3 Four wave mixing process 

Similarly, four wave mixing may also happen in different ways. One way is that in which two photon of frequency wi are 
absorbed (as pump photon) and one photon of frequency uj2 and another of frequency ws are emitted. The Hamiltonian 
representing this particular four wave mixing process is 

b'^huj2 + c+cws + g{a'^'^bc + a^h^c^). (28) 
Following the same prescription as it is used in six wave case we can write the solution as 

„2^2 

A{t) =A- 2igtA^BC + ^[AAB^BC^C - 2A^A^B^B - 2A^A^C^C - 2A'^ A^] (29) 

or 

A{t) ^ A- 2igtA^BC + gh'^[2ANBNc - NaANb - NaANc - NaA]. (30) 

The respective values of the first order antibunching (rf(l)) and second order antibunching {d{2)) can similarly be calculated 
as we have done for six wave mixing and that yields 

d(l) = {N(^Ht))-{N{t)f 

= [|a|4+52i2(_4|^|6_2|«|4)] _ [|a|4_4g2i2|^|6] (31) 

= -2g'tM'. 

and 

d(2) = (N(^Ht)) - {N{t)f 

= [|a|6-g2i2(g|^|8^g|^|6)] _ ^|^|6_6g2^2|^|8] (32) 

= -6g2i2|c,|6 

which are negative and thus they satisfy our criterion for antibunching and HOA respectively. 
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4 Second harmonic generation 

The similar procedure can be repeated for the second harmonic generation process whose Hamiltonian is 

H = hiuNi + 2hLoN2 + hg (^alaf + a\^a2^ . (33) 
The second order expression of the time evolution of annihilation operator in pump mode of second harmonic generation is 

A{t) — ai ~ 2igta\a2 + 2g^t^ ^030201 — iajaj^ . (34) 
Using the last equation along with the method used in section 2 we obtain 

= {N(^Ht))-{Nit)f 

= [|«|4+52i2(„4|^|6„2|a|4)] - [|a|4_4g2i2|^|6] (35) 



d(2) = {N(^Ht)) ~ {Nit)f 

= [|a|6-.g't'(6|a|* + 6|aH] - [|a|6-6g2i2|c,|8] (36) 

which are both negative and hence satisfies the criterion ^ for antibunching and HOA respectively. 

5 Conclusions: 

From 1I27I I32I and EEl , it is clear that all the physical systems selected for the present study show second order antibunching 
i.e. higher order sub-poissonian photon statistics. Thus the present work strongly establishes the fact that HOA is not a rare 
phenomenon. The physical systems studied in the present paper are simple and easily achievable in laboratories and thus 
it opens up the possibility of experimental observation of HOA. In case of interaction of intense electromagnetic field with 
third order nonlinear medium it was reported ^l] that the degree of antibunching {d{l)) can be tuned because they depend 
strongly on the phase of the input field which can be tuned. This is not the case with any of the physical systems studied 
in the present work. It is also clear that the higher order antibunching would not have been observed if we would have 
considered first order operator solutions (first order in g), on the other hand, if we use second order operator solutions then 
the depth of nonclassicality is found to increase monotonically with the increase of input photon number (lap). Possibly this 
monotonic increment will be ceased by the higher order perturbation terms. It is also observed that if we assume that the 
anharmonic constant and number of photon initially present in the pump mode are same for all three cases then the depth of 
nonclassicality is same in four wave mixing and second harmonic generation process and its more in six wave mixing process. 

The prescription followed in the present work is easy and straight forward and it can be used to study the possibilities 
of observing higher order antibunching in other physical systems. Thus it opens up the possibility of studying higher order 
nonclassical effects from a new perspective. This is also important from the application point of view because any probabilistic 
single photon source used for quantum cryptography has to satisfy the condition for higher order antibunching. Therefore, 
the simple prescription followed in the present work may help us to compare the existing sources of single photon. 

References 

[1] Meystre P and Sargent M III, Elements of Quantum Optics second edition Springer- Verlag, Berlin 1991. 

[2] Dodonov V V J. Opt. B. Quant, and Semiclass. Opt. 4 (2002) Rl. 

[3] Hanbury-Brown R, Twiss R Q, Nature 177 27 (1956) . 

[4] ??Hong C K and Mandel L, Phys. Rev. Lett. 54 (1985) 323. 

[5] Hong C K and Mandel L, Phys. Rev. A 32 (1985) 974. 

[6] Lee C T, Phys. Rev. A 41 1721 (1990). 

[7] Lee C T, Phys. Rev. A 41 1569 (1990). 



and 



5 



[8] M. Hillery, Phys. Rev. A 36 3796 (1987). 

[9] Giri D K and Gupta P S Optics Communication 221 (2003) 135. 
[10] An N B, J. Opt. B: Quantum Semiclass. Opt. 4 (2002) 222-227. 
[11] Pathak A and Garcia M, Communicated 



6 



